function d = fwd_TDIP(invP,etaModel)

[hs,~,~] = calcHs_TDIP(invP,etaModel,false);
d = projHs(invP,hs);